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We present a new algorithm to solve the conjugacy problem in Artin braid groups, which is faster 
than the one presented by Birman, Ko and Lee |^. This algorithm can be applied not only to braid 
groups, but to all Garside groups (which include finite type Artin groups and torus knot groups among 
others). 


1. INTRODUCTION 

Given a group G, the conjugacy problem in G consists on finding an algorithm which, given 
a,b G G, determines if there exists c G G such that a = c~^bc. Sometimes one also needs to 
compute c, for instance, when one tries to attack cryptosystems based on conjugacy in G (|Q, 

. . . . . 

We are mainly interested in Artin braid groups, which are defined, for n > 2, by the following 
presentation: 


aiaj = (Tjai (|? - j\ > 2) 

(1 ^ ^ 71 2 ) 

^Both authors partially supported by the European network TMR Sing. Eq. Diff. et Feuill. 
^Partially supported by SFRH/BD/2852/2000. 

^Partially supported by BFM-3207. 


Bn = 


CTi, (72, . 


5 ^n—1 
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The first conjugacy algorithm for braid groups was given by Garside [ pd]| . It was improved 
by Elrifai and Morton ||^ and, more recently, by Birman, Ko and Lee ([|j and Q). 

In all these algorithms, one of the key points is the existence of a finite set S C Bn, whose 
elements are called simple elements, verifying some suitable properties (we will be more precise 
later). One of the main disadvantages is the size of S, which is always greater than 3". 

In this paper we will show how one can avoid this problem by defining some small subsets 
of S, whose size is smaller than n — 1. Their elements will be called minimal simple elements. 
Unlike S, these sets of minimal simple elements are not unique for every group: The suitable set 
of minimal simple elements must be recomputed many times in our algorithm. Nevertheless, we 
will see that it is much faster to compute and use these very small subsets, than to use the whole 
S all the time. 

For instance, the known upper bound for the complexity of the Birman-Ko-Lee algorithm, to 
decide wether two braids a and b are conjugated in Bn, is 0{kPn3"') (where fc is a number that 
will be explained later, and I is the maximum of the word lengths of a and b). An upper bound 
for the complexity of our algorithm for Bn is 0{kPn'^). 

Let us mention that our algorithm, as well as the previous ones, also computes the element 
c S Bn such that a = c~^bc. Moreover, since our construction relies on the existence of simple 
elements and their basic properties, we can extend our results to a much larger class of groups, 
called Garside groups. They were introduced by Dehornoy and Paris |^. At the origin, these 
groups were called small Gaussian groups, but there has been a convention to call them Garside 
groups. They include, besides Artin braid groups, spherical (finite type) Artin groups, torus knot 
groups and others. 

One final remark: one important property of Garside groups is the existence of embedable 
monoids (for instance the monoid of positive braids, B^, which embeds in Bn). The conjugacy 
class of an element a in such a monoid is known to be a finite set, C“''(a). We will also show how 
to compute C~^{a), using the techniques mentioned above. 

This paper is structured as follows: In Section 2, we give a brief introduction to Garside 
monoids and groups; In Section 3, the known algorithms mentioned in this introduction are 
detailed; We introduce the minimal simple elements in Section 4, and in Section 5 we present 
our algorithms in detail; Gomplexity issues are treated in Section 6 and, finally, some effective 
computations are described in Section 7. 

2. GARSIDE MONOIDS AND GROUPS 

The results contained in this section are well known, and can be found in 0,0,0, i, 
, Q and 01 . We will dehne the Garside monoids and Garside groups, and explain some basic 
properties. 

Given a cancellative monoid M, with no invertible elements, we can define two different partial 
orders on its elements, ^ and Given a,b G M, we say that a ^ b {b a) ii there exists c G M 
such that ac = b (b = ca), and we say that a is a left (right) divisor of b. 

In this situation, we can naturally define the (left or right) least common multiple and greatest 
common divisor of two elements. Given a,b G M, we denote by a V 6 the left 1cm of a and b, if 
it exists. That is, a minimal element (with respect to ^) such that a aV b and b ^ aW b. We 
denote by a A 5 the left gcd of a and b, if it exists. That is, a maximal element (with respect to 
a), such that a Ab a and a Ab A b. 

Definition 2.1. Let M be a monoid. We say that x G M is an atom if x ^ 1 and if x = yz 
implies y = 1 or z = 1. M is said to be an atomic monoid if it is generated by its atoms and, 
moreover, for every x G M, there exists an integer Nx > 0 such that x cannot be written as a 
product of more than Nx atoms. 
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Definition 2.2. We say that a monoid M is a Gaussian monoid if it is atomic, (left and 
right) cancellative, and if every pair of elements in M admits a (left and right) Icm and a (left 
and right) gcd 

Definition 2.3. A Garside monoid is a Gaussian monoid which has a Garside element. A 
Garside element is an element A G M whose left divisors coincide with their right divisors, they 
form a finite set, and they generate M. 


Definition 2.4. The left (and right) divisors of A in a Garside monoid M are called simple 
elements. The (finite) set of simple elements is denoted by S. 

It is known that every Garside monoid admits a group of fractions. So we have: 

Definition 2.5. A group G is called a Garside group if it is the group of fractions of a 
Garside monoid. 


The main example of a Garside monoid (actually the monoid studied by Garside) is the Artin 
braid monoid on n strands, Bif. It is defined by Presentation considered as a presentation for 
a monoid. Its group of fractions is the braid group Bn, and Garside |I1 showed that Bif C i?„. 
Actually, every Garside monoid embeds into its corresponding Garside group |^. 

The classical choice of a Garside element for is the following: A = 

{aia 2 ■ ■ • Un-i) (o'i(T 2 • • ■ cr„_ 2 ) ■ • • ((71(72) (7i. It Can be defined as the positive braid (braid in Bif) 
in which any two strands cross exactly once (where, as usual, ai represents a crossing of the 
strands in positions i and i + I). It is represented in Figure ^ for n = 4. The simple elements 
in this case are the positive braids in which any two strands cross at most once. Then one has 
#(^) = n\ 
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FIG. 1 The Garside element A G B^. 


Another important example of Garside monoid is the Birman-Ko-Lee monoid 
the following presentation: 


which has 


BKLf = / ats(n >t>s>l) 


atsarq = arqttts if (t - r) (t - q) {s - r) {s - q) > Q 
atsdsr = atrats = asratr where n>t>s>r>l 


( 2 ) 


Its group of fractions is again the braid group i?„. The usual Garside element in BKLf is 
S = an,n-ian-i,n -2 ''' 02 . 1 - The advantage of this monoid with respect to Bf^ is that #(S) = Cn, 
where = nlfn+iy. ^ Catalan number. Hence, the number of simple elements is 

much smaller in this case, but it is still quite big, since C„ > 3”. Notice also that |5| = n — 1, 
while in i?+, |A| = . 

As we mentioned before, there are other examples of Garside groups, such as finite type Artin 
groups, or torus knot groups (see to find more examples of Garside groups). 
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From now on, M will denote a Garside monoid, G its group of fractions and A the corre¬ 
sponding Garside element. Since M C G, we will refer to the elements in M as the positive 
elements of G. 

From the existence of left Icm’s and gcd’s, it follows that (M, has a lattice structure, and S 
becomes a finite sublattice with minimum 1 and maximum A. See in Figure ^the Hasse diagram 
of the lattice of S in Bf, where the lines represent left divisibility (from bottom to top). The 
analogous properties are also verified by 





CO 


FIG. 2 The lattice of simple elements in Bf. 

Definition 2.6. For a G M we define LM (a) G S as the maximal simple left divisor of a, 
that is, LM (a) = A A a. We also define RM (a) as the maximal simple right divisor of a. 

Proposition 2.7 (0). For a G G, there exists a unique decomposition a = A^’oi---a/, 
called left normal form of a, where: 

1. p = max {r G Z : A“'’a G M} (hence ai ■ ■ ■ ai G M). 

2. Oi = LM (oi • ■ • ai) G S'\{A, 1}, for all i = I ,..., I- 
Symmetrically, one defines the right normal form of a G G, using RM. 

Sometimes, if we are dealing with elements in M and it does not lead to confusion, we will 
say that an element w = wi ■■■ Wt G M is in left normal form to express that Wi G for all 

i and, for some p > 0, the normal form of w is A^'iCp+i ■ ■ ■ Wt- 

Later we will use these technical results: 

Lemma 2.8 ([^, Prop. 2.1). Let G M be in left normal form, and G 

M in right normal form. For every v G M, one has LM(vwi---Wt) = LM (vwi) and 
RM {xi ■ ■ ■ Xtv) = RM {xtv). 

Lemma 2.9 ([^, Prop. 5.3). Let w = wi ■ ■ ■ Wt G M be written in right normal form. If we 
write w in any other way as a product oft simple elements, w = ui ■ ■ • Ut, then wi -< ui. 
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Lemma 2.10 (Q, 3.1). Letw = wi ■ ■ ■ Wt € M be written in right normal form, and let s G S. 
Then we can decompose Wi = w[w'f, for all i, in such a way that the right normal form of ws 
is (w'i){wiw' 2 ) ■ ■ ■ s) if it has t + 1 factors, or (W 1 W 2 ) ■ ■ ■ s) if it has t 

factors. 


Corollary 2.11. Let w = Wi ■ ■ ■ Wt € M be written in right normal form. Let s G S and 
suppose that we can write ws as a product oft simple elements, that is, wi ■ ■ -Wts = ui- ■ - Ut. 
Then wi -< ui. 


Proof. Since ws can be written as a prod uct of t simple elements, then its right normal form 
has t factors, say vi ■ ■ - Vt. By Lemma 2.10 , wi -< vi, and by Lemma vi ^ ui, so the result 
follows. ■ 


We end this section with a last property of Garside groups: There is a power of their Garside 
element which belongs to the center. For instance, in Bn the element = 5” generates the 
center of Bn. 


3. KNOWN ALGORITHMS FOR THE GONJUGACY PROBLEM. 

We present here the Elrifai-Morton algorithm for the conjugacy problem in braid groups 
which is also valid for Garside groups, as can be seen in 

It goes as follows: for every element a £ G, it computes a finite subset G®“™(a) of the 
conjugacy class of a. This set is shown to be independent of a, so it is an invariant of its conjugacy 
class. Therefore, two elements a and b are conjugated if and only if G®“’"(a) = G®“™(6). 

Let us explain the algorithm in more detail. 


3.1. Definition of C-™(a) and C®"'”(a) 

Proposition 3.1. 0, El Let a = A^ai ■ ■ ■ ai G G be in left normal form. Then the right 
normal form of a is as follows: a = Xi ■■ ■ xiA^, where I and p are the same as above. 

Definition 3.2. Let a = A^ai ■ --ai G G be in left normal form. We define the infimum, 
supremum and canonical length of a, respectively, by inf (a) = p, sup (a) = p + I, and ||a|| = I . 

Definition 3.3. Let a G G and denote by G (a) the conjugacy class of a. We define 
the summit infimum, the summit supremum and the summit length of a as, respectively, 
max {inf (x) : x G G (a)}, min (sup (x) : x G G (a)} and min{||a:|| : x G C (a)} . 


Definition 3.4. Let a G G. 


1. For every integer m, we define G-'^{a) = {v G G{a) : inf(z;) > m}. 

2. We define the summit class of a, G®^’” (a), as the subset of G(a) containing all elements 
of minimal canonical length. 


Remarks: 

1. One has G^°(a) = G(a) n M = G+(a). 

2. In 01, G®“™ (a) is called the Super Summit Set. 

Proposition 3.5. 00 For every b G G®“’”(a), the infimum, supremum and canonical 
length of b are equal, respectively, to the summit infimum, the summit supremum and the summit 
length of a. 


It is known that G-'^{a) and G®"*” (a) are finite sets, 
(a) ^ (j), then G™™ (a) C G^™ (a). 


Moreover, by Proposition 3.5 


if 
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3.2. Cycling and decycling 

Let T : G —> G be the automorphism defined by r(a) = A~^aA. The restriction of r to S' is 
a bijection r : S —> S. 

Definition 3.6. Let a = A^ai ■ ■ ■ ai € G be written in left normal form. The functions 
cycling and decycling are the maps c and d, from G to itself, defined by: 

c { a ) = A^a2 ■ • • a/T“^ (oi); 
d(a) = {ai)ai---ai-i. 

Notice that c(a) and d (a) are conjugates of a. Furthermore, for every a G G, inf(a) < 
inf(c(a)) and sup(a) > sup(d(a)). 

Suppose that we have an element a G G, such that inf (a) is not equal to the summit infimum 
of a. Then we can try to increase the infimum by repeated cycling. By (and p^), this 
always works: there exists a positive integer k such that inf(c^(a)) > inf (a). We know a bound 
for this integer k only for some special Garside monoids and groups: If M is homogeneous, i.e. 
it has only homogeneous relations (for instance, if M is 5+ or BKLf^), then every two words 
representing an element a G M have the same length, denoted |a|. It is shown in [p| that, in this 
case, k < |A|. 

Therefore, by repeated cycling, we can conjugate a to another element a of maximal infimum. 
Even if M is not homogeneous, we know that we reached the summit infimum when we enter 
into a loop: at some point c^{v) = v for some v conjugated to a. This always happens since the 
set G-'^{a) is finite for every m, in particular for the summit infimum. 

Once a is obtained, we can try to decrease its supremum by repeated decycling. By 
(and p^ ), this also works: either we enter into a loop, and then the supremum is minimal, or 
there exists an integer k such that sup(d^(a)) < sup(a). Again by [^, fc < |A| in homogeneous 
monoids. 

Therefore, using repeated cycling and decycling a finite number of times, one obtains an 
element a G G®"™ (a). And, if M is homogeneous, this can be done in polynomial time in |a|. 

3.3. The Elrifai-Morton algorithm 

Once that we obtained an element a G G'*“™(a), we can construct the whole G®"™ (a), by 
using the next result: 

Proposition 3.7. j^, For u,v conjugate elements in G^^'^{a) (resp. G-"^{a)), there 
exists a sequence u = ui,U 2 , ■■■,Uk = v of elements in G®“"*(a) (resp. G-'^{a)) such that, for 
1 = 1,... ,k — 1, Ui and Ui+i are conjugated by an element in S. 

The Elrifai-Morton algorithm does the following: Given a,b G G it computes, using cyclings 
and decyclings, a G G®“'"(a) and b G G®“'"(6). Then it defines Vi = {a} and it computes, by 
recurrence, 

E, = {s-^s; sGS,vG E_i}nG®“’"(a). 

Since I G S', this creates an ascending chain of subsets of G®“'"(a). By the above proposition, 
one has Vk = I4+i for some k, and then Vk = G®“'"(a). Hence, when the chain stabilises, the 
whole G®“'"(a) has been computed. Then a and b are conjugated if and only if & G (jsum 

Remark 3.8. This algorithm can be modified to compute G-"^{a) for a G M and m G Z. We 
just need to replace G®“'"(a) by G-'^{a) in the above discussion. 
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Notice that C®“"*(a) (resp. (7-™(a)) is computed at the cost of conjugating every element 
in C®“"*(a) (resp. C'-'"(a)) by every element in S. All these sets are quite big, and this makes 
the algorithm to be slow. In what follows, we will get rid of the problem caused by the size of S, 
using the minimal simple elements. 

4. MINIMAL SIMPLE ELEMENTS 

In this section we shall dehne some very small subsets of S, which will enable us to compute 
C-'^{a) and (7®“™(a), for a G G, much faster than the previous algorithms. 

Recall the definition of the partial order -< in M. 

Definition 4.1. Let V he a property for simple elements. We denote by S-p the set of simple 
elements satisfying V. The set of minimal simple elements for V, min(S''p), is the set of minimal 
elements (with respect to <) in Sp. 

We shall enforce V to be closed under g.c.d, that is, if si, S 2 S Sp then si f\ S 2 & Sp. Let us 
see that, under this assumption, the set min(iS'-p) turns to be very small. For every atom x £ M, 
let mult{x) = {s £ S'; x ^ s}. 

Lemma 4.2. Suppose that V is closed under gcd, and let x be an atom of M. If the set 
Sp C\mult{x) is non-empty, then it has a unique minimal element, that we denote px. 

Proof. Suppose that there are two distinct minimal elements si,S 2 £ Sp H mult{x). Since 
Si,S2 S Sp, then si A S2 £ Sp. Moreover, since x divides si and S2, it also divides si A S2. 
Therefore si A S2 S Sp fl mult(x), so si and S2 cannot be both minimal. ■ 

Corollary 4.3. Suppose that M hasm atoms. IfV is closed under gcd, then ff[\mTi{Sp)) < 
m. 


Proof. Notice that every element in M must be divisible by an atom. Take s £ min(Sp) and 
consider an atom x < s. Since s is minimal in Sp, it is also minimal in Sp H mult(x). Hence 
s = Px. Therefore 

min(S-p) C {px ■ x is an atom} 

and the result follows. ■ 

Example 4.4. In Bf[ there are n — 1 atoms, namely cti,... ,an-i. Therefore, if P is a 
property closed under gcd, then min(S-p) has at most n — 1 elements, while ff{S) = n! 

Example 4.5. In BKL'f there are atoms (the generators in Presentation^. Hence, 

if V is a property closed under gcd, then ff{m.m.{Sp)) < ; while ff{S) = Cn > 3". 

We must now define some suitable properties, closed under gcd, that will allow us to compute 
(7-™(a) and (7®“™(a), for a £ G. These properties will depend on some given elements in M, so 
we will have an infinite number of properties, each one corresponding to a set of minimal simple 
elements. 


4.1. Minimal simple elements to compute C7™(a) 

Definition 4.6. Let a G G and v £ G-™(a), for some m £ Z. We will say that a simple 
element s satisfies the property if it conjugates v to an element in C-"^(a), that is, s~^vs £ 
G^™(a). 
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Proposition 4.7. (Caracterization of elements satisfying V^^). If v G C-'^(a), one can 
write V = A'^w, where w G M. Then a simple element s satisfies the property if and only 

ifT™‘{s) -< ws. 

Proof. The first assertion comes from the definition of infimum. Let then v = A'^w, where 
w G M, and let s G S. One has s~^vs = s~^A^ws = A'^t'^{s~^)ws = A"*(t™(s))“^ws. Hence, 
s satisfies if and only if G M, that is, t™(s) ^ ws. ■ 

Proposition 4.8. For every v G M and every m G Z, the property is closed under gcd. 

Proof. Suppose that si and S 2 satisfy and let s = si A 32- Notice that t preserves gcd’s, 
since it preserves left divisibility. Hence r(s) = r(si) At(s 2 ), and thus t"^(s) = At"‘(s 2 ). 

One has t’^(s) -< r'"(si) ^ usi and t^(s) -< t’”(s 2 ) 'VS 2 . But it is easy to show that, for 

every v G M, vsi A VS 2 = vs. Hence, since t'^{s) divides vsi and vs 2 then it divides its gcd, i.e. 
r'"(s) A vs. Therefore, s satisfies and the result follows. ■ 


Definition 4.9. For every v G C-"‘(a), we define = min{S.p>m). That is, is the 
set of minimal simple elements (with respect to among those who conjugate v to an element 
in C^™(a). 


Notice that, by Corollary O and Proposition the cardinal of for every v G C-’^(a) 
is no bigger than the number of atoms in M. Moreover, we have the following result, analogous 
to Proposition 


Proposition 4.10. Given u,v G C'-"*(a) for some a G G, there exists a sequence u = 
ui,U 2 ,..., Uk = V of elements in C-'^(a) such that, for i = 1 , ..., fc — 1, the elements Ui and Ui+i 
are conjugated by an element in 5'^™. 


Proof. Just notice that any left or right divisor of a simple element is also a simple element, 
and then decompose every simple element in the sequence given by Proposition 3.7 into a product 
of minimal ones. ■ 


This result implies that, in order to compute C-™(a) for a G M, it suffices to conjugate every 
V G C-'^{a) by the elements in the small set 5'^'". 


4.2. Minimal simple elements to compute C®“™(a) 

Definition 4.11. Let a G G, and let v G (7®“™(a). We will say that a simple element s 
satisfies the property P®“™ if it conjugates v to an element in (7®“™ (a). In other words, if the 
canonical length of s~^vs is equal to the canonical length of v (which is the summit length of a). 


Proposition 4.12. For every v G C'®“™(a), the property P®“"* is closed under gcd. 

Proof. Let si and S2 be two simple elements satisfying P®“™, and denote s = si A S2- Write 
Si = sri for i = 1, 2, thus ri A r2 = 1. 

Suppose that inf(u) = p and ||r;|| = t. Then v = A^v', where v' G M and we can write 
v' as a product of t simple elements (but not less). Since si satisfies P®“™, one has = 

v'si = A'P TP{sf^) v'si = A'P {tP{si))~^v'si, where G M and we can write 

it as a product of t simple elements, say xi ■ ■ ■ Xt. The same happens for (t^(s2))~^t's2 G M. 

Now consider s~^vs. By Proposition it belongs to C-P{a), that is, {tP{s))~^v's G M. We 
must show that we can write this element as a product of t simple elements. Suppose this is not 
true, and write (tP(s))~^v' s = Zi ■ ■ ■ Zt+i in right normal form (it has no more than t + 1 factors 






since it is a right divisor of v's which has t + 1 factors). One has xi- ■ - Xt = (t^(si)) ^v'si = 


'{fi)) ^(tP(s)) 




Zi _ 

zi TP(r 2 )- Therefore zi r^’(ri) A rP(r 2 ) = T^’(ri A r 2 ) = ^^(l) = 1. A contradiction. 


' ■ Zt+i is in right normal form. Then by Corollary 2.11 


zt+iri. Hen ce, Zi ■ ■ ■ Zt+ixi = TP(ri)xi ■ ■ ■ Xt, and 
zi A T^’(ri). In the same way, 


Definition 4 . 13 . For every v S C®“"*(a), we define S'™™ = That is, S™™ 

is the set of minimal simple elements (with respect to among those who conjugate v to an 
element in (7®“™(a). 


As before, by Corollary 4.3 and Proposition 4.12 , the cardinal of S™™ for every t G M is no 
bigger than the number of atoms in M. Furthermore, we can adjust the algorithm by Elrifai- 
Morton to these new sets, since we have the following result, analogous to Propositions 3.7 and 
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Proposition 4 . 14 . For u,v conjugate elements in ( 7 ®“™ (a), there exists a sequence u = 
ui,...,Uk = V of elements in ( 7 ®“™ (a) such that, for t = 1 , ..., k—1, the elements Ui and rti+i are 
conjugated by an element in S®“™. 


The proof of this result parallels that of Proposition 4.10. It implies that, in order to compute 
(jsum a £ h suffices to conjugate every v G (7®“™ (a) by the elements in S®“™. 

We have then described small subsets of S which suffice to compute C-"^{a) and (7®“™(a). 
But we stili need to show how to compute these subsets. This is what we do in the next section. 


5. ALGORITHMS FOR THE CONJUGACY PROBLEM 

We shall explain in this section our algorithms to compute (7-™ (a) and (7®“™(a), given a € G. 
Let us first explain a technical algorithm, which we did not find in the literature. Let s G S' 
and V G M. We will show how to compute their 1cm s V v. More precisely, our algorithm will 
compute a simple element s' such that s V n = vs'. We must indicate that it is well known how to 
compute the 1cm and the gcd of two simple elements, as well as the normal forms of any element 
in G. 

Algorithm 1 (for computing s' such that s V u = vs'). 

1. Compute the normal form oi v = Vi ■ ■ ■ Vt- 

2. So = s. 

3. For every i = 1,... ,t, compute Si_i V Vi, and write it ViSi. 

4. Return st- 


Proposition 5.1. Let s G S and v G M. Let st be the simple element computed by Algo¬ 
rithm 1. Then s V v = vst- 

Proof. We proceed by induction on t = sup(u). If t = 1 the result is trivial, so suppose that 
t > 1 and the result is true for t — 1. Denote v' = vi ■ ■ -vt-i. We have s V v' = v'st-i, that 
is, st-i is the smallest element such that v'st-i is divisible by s. Therefore, an element r G M 
satisfies s ^ vr = v'{vtr) if and only if st-i A Vtr, and this is equivalent to St-i V Vt F VtVt, that 
is VtSt A Vtr hence Sj r. Therefore, Sj is the smallest element satisfying s -< vst, as we wanted 
to show. ■ 
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5.1. Computation of C-^{a) 


Let a G G and m S Z. As we saw in Section the main problem to compute C-^(a) is to 
compute for every v S C-^{a). 

Let then v G C-^{a) and let x be an atom of M. Consider the set S^>m n mult{x). It is 
always nonempty, since A satisfies (for every v) and is divisible by every atom (by definition 
of the Garside element). Then, by Lemma 4.2, this set has a unique minimal element, which we 
denote now r^. 

Recall that C {r^ : x is an atom}, so our first step consists of computing r^,, for every 
atom X G M. 

Let inf(z;) = p > m, so v = A'^w where w G M. Recall that, by Proposition 4.7, a simple 
element s satisfies if and only if r™(s) ^ ws. 


Algorithm 2 (for computing r^,, minimal element in Sj,>m n mult{x)). 

1. Compute the left normal form of ■ ■ -Wt- 

2. If p > TO then return x and stop. 

3. w = wi ■ ■ ■ Wt] s = X. 

4. Compute t'"(s). 

5. Use Algorithm 1 to compute s' such that V ws = wss'. 

6. If s' = 1 then return s and stop. 

7. s = ss'; go to Step 4. 


Proposition 5.2. Algorithm 2 gives an output, and it is r^,. 

Proof. First notice that, if we conjugate any p G G by a simple element, we can decrease the 
infimum of g by at most one. Hence, if p > to, inf(a;“^ua;) > m, so Vx = x and the algorithm 
gives the correct output. 

Now suppose p = m, we have computed w G M such that v = A'^w, and we need to find 
the smallest r^,, such that x < rx and T"'(rx) -< wXx- This is done as follows: we take a simple 
element s such that x -< s -< Tx (at the first step s = x). Then we use Algorithm 1 to compute 
s' such that r™(s) V ws = wss'. If s' = 1 then t^{s) -< ws, so s = and we obtain the correct 
output. Otherwise, notice that t"^(s) -< T™(ra,) ^ wvx, so wrx is divisible by t"^(s) and by ws. 
Hence wss' -< wTx, so ss' -< rx- 

Therefore, if s is not equal to Vx, the algorithm gives an element s' ^ 1 such that s ^ ss' ^ rx, 
and it starts again checking if ss' = r^,. This process must stop, since the number of left divisors 
of rx is finite, so the algorithm finds rx in finite time, (moreover if M has homogeneus relations, 
like R+ or BKL+, we find rx in at most |A| steps). ■ 


Algorithm 3 (for computing S'^™). 

1. List the atoms of M, say Xi,... , x^. Set R = cf. 

2. For 7 = 1,... , do the following: 

2a. Compute , using Algorithm 2. 
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2b. Compute Ji = {j '■ j € R and Xj -< } and 

Ki = {j '■ j > i and Xj -< r^,}. 

2c. If Ji = Ki = (j), then set i? = i? U {i}. 

3. Return {r^i : i € R}. 


Proposition 5.3. Algorithm 3 computes 

Proof. We know by Corollary |4.3| that is the set of minimal elements in {r^. : i = 
1,... , p}. We want to find a set i? C {1,.. . , v } such that = {r^. : i G i?}. Since we could 
have r^i = r^^ for some i ^ j, and we want R to be as small as possible, we define it in the 
following way: i G R ii and only if r^^ is minimal and there is no j > i such that r^. = . 

Suppose that, for some i, we already computed the elements in {!,... ,i — l}ni? (for i = 1, 
this is the empty set). Then we compute r^. and the sets Ji and Ki defined in the algorithm. If 
i ^ R^ we have two possibilities: either there is some j < i such that is a proper divisor of 
Ta;., or there is some j > i such that r^^ -< . In the first case Ji ^ 4>, and in the latter Ki ^ cf. 

Therefore, if both sets are empty, i G R. 

Using this procedure for z = 1,... , p, the algorithm computes R, thus S'^™. ■ 

Remark 5.4. Since is just the subset of {xxi ■ z = 1,... ,p} formed by its minimal 
elements, we could have computed just by comparing the r^^ ’s and keeping the minimal 
ones. We prefer to use Algorithm 3 since it is much faster to see if an atom divides an element, 
than to compare two elements, even if these two elements are simple ones. 


Remark 5.5. Algorithm 3 can still be improved in two different ways. First, we do not need 
to compute all ; if during the computation of r^^ (using Algorithm 2), we see that Xj -< s, for 
some j G R or some j > i, we can stop Algorithm 2 and increase the index i in Algorithm 3. 
Also, we do not need to compute the whole sets Ji and Ki: if we find some element belonging 
to one of them, we can directly increase the index i. We presented Algorithm 3 as above for the 
clarity of the exposition, and because these two improvements do not really change the complexity. 


Finally, let a S G and m Gh. The following algorithm works after Proposition 4.10. 


Algorithm 4 (for computing C-"^{a)) 

1. Compute the left normal form of a. 

2. Apply repeated cycling to a, to obtain a G C-'^{a) (if it exists). 

3. If d is not obtained, return (f and stop. 

4. Set V =d, V = {a} and W = (f>. 

5. Compute using Algorithm 2. 

6. For every r G do the following: 


6a. Set w = r~^vr G G-'"(a). 

6b. Compute the left normal form of w. 
6c. If zc ^ U, set U = U U {zc}. 
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7. W = WLl{v}. 

8. li V = W then return V and stop. 

9. Take a new v S V\W] go to Step 5. 


Remark 5.6. If we take m = 0, then C-^(a) = C~^(a) = C(a) H M, so this algorithm can he 
used to compute all positive elements conjugated to a. 


5.2. Computation of C®“'"(a) 

Let now a G M and v G C'®“’"(a). Let inf(u) = p and ||u|| = t (so sup(n) = p + t), and consider 
an atom x of M. Similarly to the previous case, the key point to compute C"’“'"(a) consists on 
finding the minimal element G S-psum. n mult{x) (it exists since A belongs to this set, and it 
is unique by Lemma 4.2). 


Algorithm 5 (for computing p^, minimal in Speum nmult{x)) 

1. Using Algorithm 2 compute r^,, minimal in S^>p H mult{x). 

2. s = r^. 

3. Compute the right normal form of s~^vs, say wi - ■ ■ Wk^^- 

4. If fc = t, return s and stop. 

5. s = swi; go to step 3. 


Remark 5.7. Algorithm 5 can be explained in a very natural way: first compute rx and 
the right normal form r~^vrx = wi ■ ■ -Wkl^^- If k is not minimal, that is, if k = t + 1, start 
decycling this right normal form, that is, compute the right normal form of W 2 - ■ ■ (and 

multiply rx hy wi). If this new word has not minimal canonical length, decycle again, and so on. 
The proposition below shows that this works. 


Proposition 5.8. Algorithm 5 computes px, the minimal element in Sps^m nmult{x). 

Proof. We can assume, by Proposition that we already know rx, the minimal element in 
Sj,>p n multfx). We also know that ^ px by minimality of Vx, since pf^vpx G C"’“"*(a) C 
C'-P(a). So we have an element s G S such that rx s px, s~^bs G C-^{a) and ||s“^us|| < t + 1 
(at the first step, s = rx). We compute the right normal form of s~^vs, say wi - ■ ■ Wk^^- 

There are two possible cases: either k = toik = t + l. li k = t, then ||s“^us|| = t, so s = px 
and Step 4 gives the correct output. Otherwise, s px and there exists a non trivial element 
s' G S such that ss' = Px- In this case 


Px ^vpx 


= (s') 


~^s ^vss' = (s') 


Wl ■ 


■wt+iAPs' = Ml • • -utAP, 


where ui • • • Mt A^ is in right normal form. Hence rci • ■ • Wt+i A^’s' = 


s'ui 


■ utAP ^ 


Wi 


■Wt+lT P{s') = s'ui 


■ ■ Ut- But wi • • ■ Wt+i is in right normal form, so by Corollary 2.11 


rci -< s'. Therefore rx -< swi -< px, and s is a proper divisor of swi. Also, (swi)“^u(swi) = 
W 2 ■ • ■ Wt+iA'^wi G C-P{a) and |l(swi)“^u(sr(;i)|| < f+1. So we can set s = swi, and start again. 

This procedure must stop, finding s = px, since the number of left divisors of px is finite. (In 
homogeneous monoids, the number of steps is bounded by |A|). ■ 


Algorithm 6 (for computing 
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1. List the atoms of M, say xi, ... , Xv Set R = 4>- 

2. For i = . , 1 ^, do the following: 

2a. Compute px^, using Algorithm 5. 

2b. Compute Ji = {j '■ j € R and Xj -< pxi} and 
Ki = {j ■ j > i and xj -< pxj. 

2c. li Ji = Ki = (j) then R = RU {i}. 

3. Return {px^ : i € R}. 

Proposition 5.9. Algorithm 6 computes 5'®“’". 

Proof. This algorithm parallels Algorithm 3, and it works in the same way, since S'®”™ is the 
set of minimal elements in {px^ '■ i = l,...,n}. ■ 

Finally, the following algorithm is analogous to Algorithm 4. Using the previous algorithms, 
it will then compute C'®“’”(a) for any given a € G. 

Algorithm 7 (for computing C®“'"(a)) 

1. Compute the left normal form of a. 

2. Using cyclings and decyclings, compute a € C'®“’"(a). 

3. V = a; V = {a}; W = 

4. Compute S'®"™, using Algorithm 6. 

5. For every r € S'™™, do the following: 

5a. w = r~^vr € M. 

5b. Compute the left normal form of w. 

5c. liw^V then U = U U {w}. 

6. W = W\J{v}. 

7. If U = VF then return V and stop. 

8. Take a new v G U\1U; go to Step 4. 

6. COMPLEXITY 

In this section we shall study the complexity of our algorithms, applied to several examples 
of Garside monoids, such as R+, BKL^ and Artin monoids. We do not discuss here the general 
case, since the complexity strongly depends on the way of computing normal forms, LM(a), aV6, 
a, etc. in each particular case. 

We shall give theoretical upper bounds for this complexity. In the next section we will also 
compare, with many examples in the running time of our algorithm with the one by Elrifai 
and Morton, to show that our improvement is significant in practice. 

We know that the theoretical results in this section can be improved: we are just interested 
on showing that our algorithms are much faster than the preceding ones, while sharper bounds 
for the complexity would require a deeper study. 
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6.1. The Artin braid monoid B+. 

Complexity of computing C~™(a) 

Recall the definition of B^, and notice that the relations are homogeneous, so every two 
conjugate elements in 5+ have the same word length. Let then v G B+ be of word length 1. Its 
left normal form vi ■ ■ -Vt verifies t < 1. 

Algorithm 1. It computes the normal form of u, which takes time 0{Pnlogn) (see[M). Then 
it computes t words, si,... , St- The computation of each one takes time O(nlogn) (pq), so all 
of them are computed in time 0{ln\ogn), and the whole algorithm has complexity Ojlrnlogn). 

Now suppose that v G Bn, and its left normal form is A^wi ■ ■ -wt- One has t < 1. In Bn, the 
homomorphism r can be defined as follows: r(cri) = an-i- Hence = id, so for every word w, 
T™'{w) can be computed in time 0(|w|). 

Algorithm 2. First it computes the left normal form of v in time 0{Pn log n). li p = m, it runs 
the loop consisting of Steps 4-7. The two important steps are the following: 

• Step 4: Notice that, every time the loop is repeated, we already know t’"(s) for the 

old value of s. So in order to compute = T™(s)r'"(s'), we just need to compute 

which is 0(|s'l). Since the product of all possible s' is still a simple element, Tx, all 
repetitions of this step can be made in time 0{n{n — l)/2), which is the length of A. 

• Step 5: Here, when we apply Algorithm 1, we compute the normal form of ws, and the 
elements si,... ,St, where s runs over an ascending chain of divisors of A. As above, 
all these computations together require the same number of operations as just applying 
Algorithm 1 to wTx (see [Q). Moreover, we have computed the normal form of w at the 
beginning of Algorithm 2 so, again by fl^ , the normal form of wTx can be computed in 
time 0{lnlogn). Hence all repetitions of this step can be made in time 0{ln\ogn). 

Therefore, the theoretical complexity of Algorithm 2 is 0{Pri?). 

Algorithm 3. The only non-negligible step is the loop of Step 2, which is repeated n — 1 times 
(the number of atoms in B^), and does the following: it computes (which takes time 0{P'n?')) 
and it verifies at most n — 2 times if an atom divides Tx^ (this takes time 0{n) by |^). Hence, 
the complexity of Algorithm 3 is 0{Pn^). 

Finally, let a G be given as a word of length I in the generators (7^. 

Algorithm 4. It starts by computing the normal form of a (time O(l^nlogn)). Then it finds 
a, which takes 0{Pn^) by Q]. Next it starts a loop, which is repeated k times (the number of 
elements in C'-"*(a)), and does the following: First, it computes S'^™, taking time 0{Pn^). Then 
it runs another loop, repeated at most n — 1 times, which works at follows: 

• It computes the normal form of a word of length I, thus taking time 0{Pn\ogn). 

• It verifies if an element is in a list V, taking a negligible time compared to the previous 
step. 

Finally, it verifies if F = W, but since kF C F we just have to compare the lengths. The time 
to do this is negligible. Therefore, the complexity of Algorithm 4 (i.e. the complexity of computing 
C'-™(a)), is O (l^nlogn + + k{Pn^ + {n — I)Pnlogn)), which yields the following. 

Proposition 6.1. Given a G Bn as a word of length I, the complexity of computing C-''^{a) 
(for the Artin presentation) is O {kPn^^, where k is the number of elements in C-"'{a). 
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Remark that if we try to compute C-'^{a) using the techniques of Elrifai and Morton, we 
should use Algorithm 4, but replacing by the whole S', which has cardinality nl. The time 
would be in this case O (A:Z^(n!)nlogn). 

Complexity of computing C®“™(a) 

The study of the complexity of Algorithms 5, 6 and 7, is very similar to that of Algorithms 2, 
3 and 4. 

Algorithm 5. It starts by computing r^, taking 0{Pn^). Next it computes a right normal form 
(O(Pnlogn)), and then it does a number of decyclings, which is bounded by ®^ch 

decycling takes time 0{ln), so the whole complexity of Algorithm 5 is 0{l^n^). 

Algorithm 6. It does the same as Algorithm 3, but it computes px^ instead of Vx^ in Step 2a. 
Hence its complexity is 0{Pri^). 

Algorithm 7. It has two main differences with respect to Algorithm 4. It computes a instead 
of o (but this can be made in 0{Pn^) by Q]), and it computes S'™™ instead of S^™. Therefore 
one has the following. 

Proposition 6.2. Given a G as a word of length I, the complexity of computing C'®“™(a) 
(for the Artin presentation) is 0{kP‘rA), where k is the number of elements in C'®“™(a). 

Remark that, if we compute the complexity of the algorithm by Elrifai and Morton, using the 
above methods, we obtain O (kP {ni)n log n). 

6.2. The Birman-Ko-Lee monoid BKL+. 

Complexity of computing C’*’(a) 

We just need to follow here the same reasoning that in the previous subsection, taking into 
account the differences between BKL'f and R+. The complexities of the basic computations in 
BKLfl can be found in 0. For instance, in BKL)^, the computation of each Si in Algorithm 1 
takes time 0{n), the length of the Garside element is n — 1, and the normal form of a word 
w is computed in time 0{\w\'^n). This implies that Algorithm 1 and Algorithm 2 both have 
complexity 0{Pn). 

In order to study Algorithm 3, we must know that the number of atoms in BKLfl is ^ 

and to check if an atom divides a simple element takes time 0(n), so Algorithm 3 has complexity 
0{Pn^). 

Finally, the computation of a takes time 0{Pn^) (see 0), so the above method to compute 
the complexity of Algorithm 4 yields the following: 

Proposition 6.3. Given a G BKLn as a word of length I, the complexity of computing 
C-'^{a) (for the Birman-Ko-Lee presentation) is O ikPn^'), where k is the number of elements 
in C^'^{a). 

As above, the complexity of the previously known algorithm is much worse: 0{kPCnn), where 
3” < C„ < 4". 

Complexity of computing C®“™(a) 

We just need to know that the computation of d takes time 0{Pn^) (see Q). Hence the 
complexities of Algorithms 5, 6 and 7 are respectively 0{Pn^)^ 0{Pn^) and 0{kPn^). Therefore, 
one has: 
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Proposition 6.4. Given a € BKL^ as a word of length I, the complexity of computing 
Csum^a'j (for the Birman-Ko-Lee presentation) is 0{kPn^), where k is the number of elements 
in (7®“™(a). 

Notice that the complexity of the known algorithm was 0{kPCnn), so our algorithm improves 
it considerably. 

One interesting remark is that our algorithm works faster, a priori, for the monoid Bif than 
for BKLf. This is due to a simple fact: in our algorithm the number of atoms is more relevant 
than the number of simple elements. In BKLf, the number of simple elements is much smaller 
than in Bf^, but the number of atoms is , while in Bf^ is n — 1. 

6.3. Artin monoids 

As we mentioned in the introduction, the Artin groups of finite type are Garside groups, so we 
can apply our algorithms to the corresponding Artin monoids (see for an introduction to Artin 
monoids and groups). In we can find algorithms to deal with Artin monoids: computation 
of normal forms, greatest common divisors, division algorithms, etc. Although these algorithms 
seem to be exponential in the length of the words involved, in |Q] it is shown that finite type 
Artin groups are biautomatic, so there are quadratic algorithms to compute all of the above. 

Nevertheless, since we are mainly interested in comparing our algorithms with the previous 
ones, we just need to know the length of the Garside element A, and the number of simple 
elements in any given Artin group. Let then G be an Artin group of rank n, that is, A„, i?„, 
Dn, En (if n = 6, 7,8), (if n = 4), (n = 3,4) or hip) (if n = 2), and let h be its Goxeter 
number. It is known that |A| = where h = 0(n), and that #(S') > n!. 

Hence, if the complexity of the conjugacy algorithm by Elrifai and Morton is 0(a:n!) for some 
X depending on n and I, our algorithm will have complexity 0{xn^). This is shown by using the 
same arguments as in the previous subsections. 

7. EFFECTIVE COMPUTATIONS 

7.1. Comparison with the Elrifai-Morton algorthim 

In the previous section, we found theoretical upper bounds for the complexity of our algo¬ 
rithms. We showed that our algorithm is, in theory, much better than the Elrifai-Morton one 
(for n > 5). In this section we effectively compare the two algorithms, in the following way: For 
given n and /, (3 < n < 5 and 10 < Z < 20) we took 5000 random pairs of positive braids in 
Bn of length I (using Artin presentation), we tested conjugacy using both algorithms, and we 
compared the Average Running Time (ART) and the Maximum Running Time (MRT). We did 
the same for n = 6 and I = 10, for 1144 pairs. 

We can conclude that our algorithm is faster for n > 4, and much faster for n > 5 (We were 
not able to compute the cases n = 5 and I = 19, 20 using the Elrifai-Morton algorithm since the 
computations were too long). 

In the tables below one can see the results: We wrote F-GM for our algorithm and E-M for 
the Elrifai-Morton one. The time is given in seconds. 


n = 3 


1 

10 

11 

12 

13 

14 

15 

ART F-GM 

0.1526 

0.2011 

0.2361 

0.3038 

0.3386 

0.3951 

ART E-M 

0.1144 

0.1460 

0.1692 

0.2133 

0.2367 

0.2723 

MRT F-GM 

2.429 

3.599 

4.680 

6.080 

7.450 

6.960 

MRT E-M 

1.659 

2.539 

3.220 

4.089 

5.029 

4.599 
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1 

16 

17 

18 

19 

20 

ART F-GM 

0.3896 

0.5021 

0.5473 

0.6494 

0.7292 

ART E-M 

0.2710 

0.3392 

0.3710 

0.4329 

0.4841 

MRT F-GM 

11.299 

10.530 

12.469 

15.090 

16.539 

MRT E-M 

7.219 

6.729 

7.970 

9.950 

11.039 


71 — 4 , 


1 

10 

11 

12 

13 

14 

15 

ART F-GM 

0.3559 

0.4796 

0.6772 

0.7870 

1.0264 

1.2599 

ART E-M 

0.6118 

0.7233 

1.0127 

1.2086 

1.5909 

1.9538 

MRT F-GM 

8.680 

11.390 

16.519 

23.949 

33.969 

42.029 

MRT E-M 

16.319 

22.329 

28.440 

41.579 

61.790 

74.999 


1 

16 

17 

18 

19 

20 

ART F-GM 

1.4548 

1.7436 

2.2029 

2.6616 

2.9942 

ART E-M 

2.3106 

2.7995 

3.5548 

4.3280 

4.7226 

MRT F-GM 

41.910 

62.940 

72.940 

103.470 

148.989 

MRT E-M 

70.039 

107.969 

137.720 

173.060 

245.740 


n = 5 


1 

10 

11 

12 

13 

14 

15 

ART F-GM 

1.0997 

1.8463 

2.7657 

3.7962 

3.8195 

4.4797 

ART E-M 

7.8690 

11.1207 

17.1455 

23.2491 

26.2595 

29.7934 

MRT F-GM 

21.239 

46.070 

65.530 

88.940 

139.180 

155.260 

MRT E-M 

177.489 

322.039 

456.669 

611.609 

1068.970 

1178.790 


1 

16 

17 

18 

19 

20 

ART F-GM 

5.6410 

7.1540 

8.8198 

9.4597 

10.6614 

ART E-M 

38.7974 

51.0028 

62.0018 



MRT F-GM 

254.770 

411.320 

401.409 

516.119 

532.469 

MRT E-M 

2116.239 

3221.880 

3218.93 




n = 6 


1 

10 

ART F-GM 

2.2450 

ART E-M 

506.224 

MRT F-GM 

43.935 

MRT E-M 

7495.288 


7.2. Exhaustive computation of conjugacy classes and summit classes 

In the previous section, we saw that the complexity of all our algorithms depends on the size 
of the sets C'-™(a) or C'®“’"(a), for a G M. In the cases of B+ or BKL^, the only upper bounds 
known for these sets are exponential in n and in ^ = |a|. Nevertheless, we have the following 
(recall that, in this case, C+(a) = C-^{a) = C{a) fl B^): 
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Conjecture: (Thurston, |^) Let n be a fixed integer and let a G , having word length 1. 
There is an upper bound for C^{a) which is a polynomial in 1. 

The existence of this upper bound for (7+(a), or even for would imply the following: 

Conjecture: (Birman, Ko and Lee, [Q) For every fixed integer n, there exists a solution for the 
conjugacy problem in Bn, which is polynomial in the word length of the elements involved. 

In order to have some numerical evidence to support these conjectures, we have computed, 
for n = 3,... ,8 and several values of I, all the conjugacy classes of words of length I in 5+, as 
well as the corresponding summit classes. In the tables below we present the following data, for 
the set Wi of elements in having word length /: 

• CC+: The number of Conjugacy Classes in Wi C i?+. 

• maix C+: The size of the biggest one. That is, the number of elements in the biggest 
C+(a), for a GWi- 

• max C®“™: The size of the biggest summit class. 

• v: A representative from one of those biggest summit class. That is, an element v G 
C®“’”(a), where C®“’”(a) has maximal size. 


n=3 

/ 

CC+ 

max C~^ 

maxC®“’” 

V 

4 

3 

6 

2 

crj* (72 

5 

3 

10 

6 

afa^ 

6 

5 

12 

8 

afa'i 

7 

5 

16 

10 

0-1^2 

8 

8 

20 

12 

a^a'i 

9 

9 

29 

14 

aia'i 

10 

13 

30 

16 

afa'i 

11 

16 

40 

18 


12 

27 

48 

20 


13 

33 

64 

22 

erf cr^ 

14 

50 

80 

24 

erfer^ 

15 

70 

125 

26 

erfer^ 

16 

107 

126 

28 

erfer^ 

17 

153 

160 

30 

erfer^ 

18 

241 

192 

32 

erfer^ 

19 

349 

256 

34 

ai^a'i 

20 

542 

320 

36 

erfer^ 
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n=6 

H 

CC+ 

max (7+ 


V 

4 

13 

36 

16 

(TiCT2Cr3Cr4 

5 

22 

56 

30 


6 

44 

120 

36 

erf (72 CTs 

7 

76 

272 

72 


8 

148 

412 

124 


9 

276 

576 

208 


10 

573 

1032 

372 



n=7 

1 

CC+ 

max (7+ 

max (7®“™ 

V 

4 

14 

60 

24 


5 

26 

84 

60 


6 

56 

160 

72 


7 

104 

408 

108 


8 

215 

824 

192 


9 

424 

1160 

416 


10 

914 

1992 

744 
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n=8 

1 

CC+ 

max 

maxC'™’” 

V 

4 

15 

100 

48 

0’iCr20'3(T5 

5 

29 

144 

100 

CTlCr20’i(J4 

6 

66 

216 

144 

(Ti(720’lCT3Cr5 

7 

130 

544 

168 

(Ti(J2CriCr3(J4(Tg 

8 

281 

1236 

360 

(Ti(J2Cri(T^CT^ 
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